To be presented at: 15th Biennial ASME Vibrations Conference 

Boston, Massachusetts N95- 32692 

September 17-21, 1995 

O V~ 


DYNAMIC MODELLING AND RESPONSE CHARACTERISTICS 
OF A MAGNETIC BEARING ROTOR SYSTEM 
WITH AUXILIARY BEARINGS 





April M. Free and George T. Flowers 

Department of Mechanical Engineering 
Auburn University 
Auburn, AL 36849 

Victor S. Trent 

Department of Electrical Engineering 
Auburn University 
Auburn, AL 36849 


ABSTRACT 

Auxiliary bearings are a critical feature of 
any magnetic bearing system. They protect the 
soft iron core of the magnetic bearing during 
an overload or failure. An auxiliary bearing 
typically consists of a rolling element bearing 
or bushing with a clearance gap between the 
rotor and the inner race of the support. The 
dynamics of such systems can be quite complex. 
It is desired to develop a rotordynamic model 
which describes the dynamic behavior of a 
flexible rotor system with magnetic bearings 
including auxiliary bearings. The model is 
based upon an experimental test facility. Some 
simulation studies are presented to illustrate 
the behavior of the model. In particular, 
the effects of introducing sideloading from the 
magnetic bearing when one coil fails is studied. 
These results are presented and discussed. 

NOMENCLATURE 
C = damping, N-s/m 

Cl = clearance of auxiliary bearing, m 
D — nominal gap thickness, m 
Fx = external force vector acting on the rotor 
in 

X direction, N 

F y = external force vector acting on the rotor 
in 

Y direction, N 

G = gravitational acceleration, m/s 2 


I = rotor inertia matrix 

t = current, amp 

I( — stiffness, N/m 

k = gain value 

L = equivalent circuit length 

M — mass, kg 

m = total number of nodes 

N = total number of modes considered 

Nbi = node number at rotor left end 

Nb2 = node number at AMB 

Nb 3 = node number at auxiliary bearing 

Ni m b = node number at imbalance location 

Nt = number of turns of wire per coil 

Qx = rotor modal coordinate vector in X 

direction 

Qy = rotor modal coordinate vector in Y 
direction 

R r = radius of rotor journal, m. 

Rs = radius of auxiliary bearing bore, m. 
t = time, s 
tf = fail time, s 
v = voltage, volt 

v r = relative velocity at auxiliary bearing 
rotor/stator contact point 
Xfi = rotor physical coordinate vector in 
X direction 

Xs = stator physical coordinate vector in 
X direction 

Yfl = rotor physical coordinate vector in 

Y direction 

Y 5 = stator physical coordinate vector in 

Y direction 

a = acceleration, rad/sec 2 
6 = dynamic clearance 



r = 

$ = rotor free-free modal displacement matrix 
= rotor free-free modal rotation matrix 
rp = imbalance vector 

= permeability of free space 
fi r = permeability of silicon steel 
SI = rotor operating speed, rad/s 
r = time constant 

0 = angular position of the shaft 
£ = modal damping coefficient 

Subscripts 
a = first 

B = auxiliary bearing 
b = second 
C = contact 
g — gain constant 
h = horizontal 
/ = left bearing 
r = rate constant 
v = vertical 
x = x — direction 
y = y — direction 
= bias 

1 = top magnet 

2 = right magnet 

3 = bottom magnet 

4 = left magnet 

Superscripts 

a = current amplifier 
c = controller 
l — lead 

p — sensor amplifier 
INTRODUCTION 

In recent years, the use of active magnetic 
bearings (AMB) for turbomachinery support 
has been an area of interest in both academia 
and in industry. Magnetic bearings provide 
the potential for significant improvements over 
other types of rotor supports, including elimi- 
nation of wear, bearing friction-related energy 
losses and a means of actively suppressing ro- 
tor vibration. A critical feature of any mag- 
netic bearing supported rotor system is auxil- 
iary bearings to protect the soft iron core from 
rotor contact (and subsequent damage) during 
an overload or failure of the AMB. The present 
work is concerned with developing a suitable 
model and investigating the dynamic behavior 
of such a system. 


EXPERIMENTAL TEST FACILITY 

The model and control system development 
presented in this study is for an experimental 
test facility. The apparatus consists of a radial 
magnetic bearing that is supporting the right 
end of a rotor. The left end of the rotor is 
supported by a ball bearing suspended in a 
frame by four springs. A photograph and the 
corresponding schematic diagram are shown in 
Figure 1. 

The apparatus has two basic components: 
the rotor shaft and the magnetic bearing. The 
shaft is made of steel and is 0.0098 meters in 
diameter and 0.4572 meters in length. A steel 
disk of diameter 7.62 cm, thickness 2.54 cm, and 
mass 0.810 kg is placed at the midpoint of the 
bearing span. Threaded holes on the disk allow 
for imbalance to be added to the system. The 
rotor is driven by a variable speed motor with 
a controller. Shaft vibration is measured using 
eddy current proximity displacement sensors 
fixed to measure displacement in both the 
vertical and horizontal directions. 



Figure l.a Experimental Test Facility 

The magnetic bearing is based upon a design 
described in Humphris, et al (1988). The basic 
parameters for the bearing are shown in Table 
1. The bearing consists of four electromagnets 
equally spaced around a soft-iron core. A pho- 
tograph and schematic diagram of the bearing 
are shown in Figure 2. 

SIMULATION MODEL 

A simulation model was developed lor the ro- 
tor system described above. 1 he model has 
three principal components the rotor, mag- 
netic bearing, and electronics. The governing 
equations for each are shown below. 





Figure l.b Schematic Diagram of Experimental 
Test Facility 



Figure 2. a Magnetic Bearing Assembly 



Figure 2.b Schematic Diagram of Magnetic 
Bearing 


The rotor is modelled using the free-free 
bending mode shapes and natural frequencies 
obtained through finite element analysis. The 
finite element code uses 1 1 stations and the first 
four modes (two rigid body and two flexible 
modes) are included in the simulation model. 
The rotor equations of motion can be expressed 
in terms of modal coordinates as 

Qx + ftrQ„ + u^Qx + * T F r = 0, (1.0) 


Q v - nrQx + u£q v + * t f v = o, (i.6) 

The auxiliary bearing consists of a ball bear- 
ing with a clearance between the rotor and the 
stator. The governing equations are 

MbX&-\-C bxX b 4* KbxXb = F Xtaux ,(l.c) 
^B^B+^By^B + ^By^B = Fy rOU ar> (1**0 

where 

F x ,lin(N bl ) = $ATh(#(xQx + C, x Qx) 
F v ,k„(JV h ) = *N„{K,y Q y + C,y Qv) 
Fx,imfc(Af<m6) = tffocose + a^ 6 sin0, 

F y,imb{N im b) = ft 2 i(} 6 sin0 - arj) 6 cos0, 
Fx,a«x(lVfc3) = *n„(F NCOS0 - F Rsinfl) 
Fy iat , x (JVfc 3 ) = ^N t3 (F N sin/3 + Frcos/3 ) 
Fgrav = * T MG 
F s = 4>Kc(6 - D). 

If t? r = 0, then F/i = required static friction 
force to satisfy this condition. If v r / 0 or 
if required static friction force > /ZjF^, then 

Ffl = /ifcFjsr- 

F x = F Xf /^ n + F X| tm6 F x>aux — F Xf a m 6 

Fy = Fy ( ^ n + Fy ||m 6 + Fy >aux “ F y, fl m6 + F gmV 


where 







implemented using analog circuitry. The con- 
troller is modelled as 


sin(/?) = 


Yr-Ys 

v/(Xr-Xs) 2 + (Yr-Y 5 ) 2 , 


cos((3) = 


Xr-Xs 

^(Xr-X s ) 2 + (Yr-Y s ) 2 , 


S=y/(* Qx-X t )2 + ($Q y -n) 2 , 

<t> = 1 if 6 > A 
0 otherwise 
Qx = * T X r , 

Qv = * T Y r , 


with 


x r = {X rU X r2 , . . . , X rm } T , 

Y r = {Y rl ,Yr2,...,Y Tm } T . 


dr, 

dt 


? = k a v>r c + (k g + k r )f t -T*v e v ( 3.a) 


dr ( c 

dt 


K -k g vlr c + (k g + k r )f t -r c v c h (3.6) 


The control circuitry also includes a lead 
network. The governing equations are 


dv[, 

dt 


dt 




(4.a) 


dv‘ h 1 c dv c h 1 , ... 

-i = rf”* + IT - 7i Vh (4 t) 


The physical displacements of the rotor at 
the bearing and imbalance locations can be ob- 
tained using the following coordinate transfor- 
mation: 

N 

X r k — ^ ^ *ki Qxi't 

t=l 

N 

Yrk = 

1=1 

(k = , Nb2) N 63, N tmb) 

The position of the shaft at the magnetic 
bearing location is measured using horizontal 
and vertical proximity sensors. They are mod- 
elled as linear first order systems. 


dvl __ k v y - v v v 
dt r p 


(2.a) 


Four current amplifiers (one for each coil) 
supply current to the magnetic bearing. They 
are modelled as 


= . J_ i? + ±*;„i 

dt r° 1 r° 1 


(5.a) 


dt 


1 . 1 


^ = -—*5 + — *?t»i 

dt r° 2 r° 2 v 


(5.6) 


di | 
dt 


^ .a . ^ l«..I 

7^*3 + 


(5.C) 


dt t* 4 r a 4 h 


(5.d) 


du p Ar p x - v p 
dt r v 


( 2 . 6 ) 


The active control system is based upon a 
proportional-derivative (PD) control law and 


The forces supplied by the magnetic bearing 
to the rotor are assumed to be decoupled in the 
horizontal and vertical directions. The effects of 
flux saturation are included in the model based 
upon the approach described by Lewis (1993). 
The equations for these forces are 



AiV'j' ( — Y, T 

( 2 (D- x „ B3 )+-) 2 

hr 

— ) 

(2(D + x Nb7 )+—) 2 

rr 

(6.a) 




*1 

(2 (D - DNb?) + — f 

hr 

^1 _) 

(2 (D + Vn B2 ) + ~Y 


Mr 


( 6 . 6 ) 



Figure 3 Block Diagram of Simulation Model 


DISCUSSION 

Using the rotordynamic model described in 
the preceding paragraphs, a study of the com- 
bined dynamics of a magnetic bearing sup- 
ported rotor with auxiliary bearings was per- 
formed. The parameters used in this work were 
identified from experimental evaluation of the 
test rig and are shown in Table 1. All of the 


Parameter 

Value 

Units 

K| x , K| y 

17,510 

N/m 

Clxj C|y 

2.0 

N-s/m 

K c 

87,557 

N/m 

Cc 

0 

N-s/m 

Kb 

17,510 

N/m 

Cbx» Cj3y 

2.0 

N-s/m 

M b 

0.033 

kg 

A 

3.0 x 10-“ 

kg-m 

h? 

7.8 x 10 3 

“ 

kg 

13 

- 


23 

- 

k a 

0.45 

amp/v 

rP 

1.59 x 10 -b 

- 

rT 

1.8 x 10“^ 1 

- 

Tfc 

1.59 x 10“ 4 

- 

r° 

5.31 x 10~ & 

- 

A 

3.42 x 10" 4 

m 2 

N t 

164 

- 

D 

0.9 x 10 -3 

m 

ho 

1.26 x 10"“ 

- 


4.2 

volt 

< 

1.9 

volt 


1.0 

volt 

5 ? 

1.9 

volt 


Table 1 Simulation Model Parameters 


responses illustrated in the figures are for this 
base parametric configuration, unless otherwise 
indicated. Figure 3 shows the variation of sys- 
tem natural frequencies with rotor speed for the 
base parametric configuration (without auxil- 
iary bearing contact). The rotor speed of 200 
rad/sec was selected as a reasonable value near 
a critical speed in order to illustrate a worst 
case behavior of the system. 

Figure 4 shows the rotor response under nor- 
mal operating conditions. The horizontal and 
vertical responses are of the same amplitude 
and frequency at steady state. The response 
consists of a transient region of about 0.5 sec- 
onds in length followed by steady synchronous 
oscillations. The initial conditions are set so 
that the rotor contacts the auxiliary bearing 
during the transient phase of the oscillations. 
The rotor vibrates downward, contacts the aux- 
iliary bearing, bounces several times, and then 
settles into a steady synchronous whirl of about 






















0.05 mm in amplitude. 

Of particular interest in this study is the use 
of the remaining coils of the AMB (when one 
fails) to provide a sideloading force to reduce 
rotor vibration. The test configuration consists 
of allowing coil 2 to fail and investigating the ef- 
fects of siueloading on the system. Coils 1 and 3 
are kept in the same electronic configuration as 
if no failure had occurred. However, the voltage 
in coil 4 is set to a value that will pull the shaft 
horizontally over to the auxiliary bearing. Two 
basic scenarios for rotor control during failure 
are considered. Figures 5. a - 5.c illustrate the 
rotor behavior under the conditions for the first 
scenario, when the rotor speed is maintained at 
a constant value. For each case, the AMB fails 
at tj= 0.5 seconds. There is some perturbation 
of the vertical motion of the rotor when the hor- 
izontal coil fails due to coupling of the horizon- 
tal and vertical motion (primarily) through the 
gyroscopic terms. However, the steady state 
vertical responses are not significantly affected. 
Figures 6.a - 6.c show the rotor behavior for 
the same failure conditions (with sideloading) 
as above for the second scenario, when the ro- 
tor is decelerated at a constant rate upon AMB 
failure. Again the response amplitudes tend to 
be lowest for smaller clearance values. It is ob- 
served that the responses may be quite high as 
the rotor contacts the auxiliary bearing, even 
with sideloading from the remaining coil. How- 
ever, sideloading is shown to be quite effective 
in encouraging rotor contact that will dissipate 
the vibrational energy. The response ampli- 
tudes tend to decrease as the auxiliary bearing 
clearance is decreased. However from a practi- 
cal perspective, there are definite design trade- 
offs that must be considered. If the clearance 
is too small, the rotor may strike the auxiliary 
baring during normal operation of the mag- 
netic bearing if the rotor is perturbed. This 
tends to unnecessarily shorten the life of the 
auxiliary bearing. If the clearance is too large, 
the rotor vibration (upon AMB failure) will be 
excessively high ana may damage the rotor or 
the magnetic bearing structure. 

The effect the sideloading voltage on the sys- 
tem dynamics is particularly interesting. The 
magnetic bearing forces are nonlinear functions 
of the dynamic clearance and the current. If one 
coil fails and the remaining coils are used to pro- 
vide sideloading, the natural approach might 
seem to issue a constant command voltage, with 
the objective of having a constant sideloading 
force at steady state. However, the inherently 
unstable behavior of the magnetic bearing with- 
out active control coupled with the vibrational 



rotor speed (rad/sec) 


Figure 3 Rotor Natural Frequencies as 
a Function of Rotor Speed 
(for the parameters of Table 1 ) 
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Figure 4 Rotor Response Without AMB Failure 
ft = 200rad/sec, Cl, = ^ 


effects between the auxiliary bearing and the 
rotor produces a severely unstable response. A 
remedy is to command a voltage that is pro- 
portional to the horizontal shaft displacement. 
That is, v a = v a a {D + x N b? ), where v a Q is a con- 
stant. With an ideal amplifier, any value of v 0 
can be used with such a system. However due 
to the delay in the dynamics from the time con- 
stant of the power amplifier, unstable behavior 
may also result if v a is too large. Figure 7 shows 
the maximum that can be commanded with 
stable behavior as a function of power amplifier 
time constant. These results were developed by 
fitting a smooth curve through data points ob- 
tained by observing the behavior of the simula- 
tion model for various sideloading voltages and 
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|6.a Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = ^p 
t j = 0.5, a = 0 



0.5 1.0 

Time (sec) 


i.b Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = y 
ty = 0.5, a = 0 



Figure 5.c Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = y 
tj - 0.5, a = 0 



Time (sec) 


Figure 0.a Rotor Response With AMB Failure 
ft = 200rad/sec, Cl = ^p 
t f = 0.5,a = -100 


I time constants. It must be concluded 
pr amplifier dynamics are a critical fac- 
Tselection of appropriate levels of bear- 
ing sideloading. 


CONCLUSION 
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A simulation model has been developed for a 
magnetic bearing supported flexible rotor sys- 
tem with auxiliary bearings including frictional 
effects. The model has been described in detail 
with experimentally obtained model parame- 
ters. The response characteristics for a vari- 
ety of system configurations were presented and 
discussed, including the effects of sideloading. 
Some guidelines were given for the selection of 
appropriate levels of sideloading. 


REFERENCES 

Gondhalekar, V., and Holmes, R., “Design 
of a Radial Electromagnetic Bearing for the 
Vibration Control of a Supercritical Shaft,” 
Proceedings of the Institution of Mechanical 
Engineers, Vol. 198C, No. 16, pp. 235-242. 

Habermann, H., and Liard, G., “An Active 
Magnetic Bearing System,” Tribology Inter- 








Figure 6.b Rotor Response With AMB Failure 
SI = 200rad/sec, Cl = y 
tf = 0.5, a = —100 



V (sec 1 ) 


Figure 7 Maximum Sideloading Force as 
a Function of Power Amplifier 
Time Constant, r a 
(for parameters of Table 1) 



Figure 6.c Rotor Response With AMB Failure 
fl = 200rad/sec,C£ / = y 
tf = 0.5, a = -100 


E., “Magnetic Bearing Control Systems and 
Adaptive Forced Balancing,” IEEE Control 
Systems , Vol. 14, No. 2, pp. 4-13. 

APPENDIX LINEAR CURRENT AND 
POSITION STIFFNESSES 

For design purposes, a linearized model for 
the magnetic bearing forces is needed. 

The current stiffnesses are: 
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1 2D 2 
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■i _ ppAk^v 2*Mt 
2 “ 2D 2 
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3 2D 2 


(7.c) 


PoAk*v b 2 'Nl 

A < ~ 2D 2 

The position stiffnesses are: 



2D 3 


(7.d) 


(8.n) 


K 


P 

2 


2D 3 


( 8 . 6 ) 






